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limits of validity of these approximations. For mu approximately less than 
10 iieV/G (E approximately less than 10 keV at L equivalent to 3) the ion 
drift period on the final < r ing-current; ) drift shell of interest (L 
equivalent to 3> exceeds the duration of the main phase of our model 
storm., and we find that the transport of ions to this drift shell is 
appropr lately idealised as direct convective access, typically from open 
drift paths. Ion transport to a final closed drift path from an open 
(plasma-sheet) drift trajectory is possible for those portions of that 
drift path that lie outside the mean stormtime separatrix between closed 
c\nd open drift trajectories, For mu approximately 10-25 MeV/G (110 keV 
approximately less than E approximately less than 280 keV at L equivalent 
to 3) the drift period at L equivalent to 3 is comparable to the 
postulated 3-hr duration of the storm, and the mode of transport is 
transitional between direct convective access and transport that resembles 
radial diffusion.. (This particle population is transitional between the 
ring current and radiation belt)., For mu approximately greater than 25 
MeV/G (radiation-belt ions having E approximately greater than 880 keV at 
L. equivalent to 3) the ion drift period is considerably shorter than the 
main phase of a typical storm, and ions gain access to the ring-current 
region essentially via radial diffusion. By computing the mean and 
mean-square cumulative changes in 1/L among (in this case) 18 
representative ions equally spaced in drift time around the steady-state 
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drift shell of interest (L equivalent to 3). we have estimated (from both 
our forward and our time-reversed simulations) the time-integrated 
radial-diffusion coefficients D(sup simMsub LL) for particles having 
selected values of mu approximately greater than 15 MeV/G. The results 
agree surprisingly well with the predictions (D(sup qlHsub LL) ) of 
quasilinear radial diffusion theory-, despite the rather brief duration 
(approximately 3 hrs) of our model storm and despite the extreme 
variability (with frequency) of the spectral-density function that 
characterizes the applied electric field during our model storm. As 
expected., the values of D(sup simXsub LL.) deduced ( respect i vely ) from our 
forward and time-reversed simulations agree even better with each other 
and with D ( sup simXsub LL) when the impulse amplitudes which characterize 
the individual substorms of our model storm are systematically reduced,. 
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Stormtime Transport of Ring Current and Radiation Belt Ions 

Margaret W. Chen, Michael Schulz, Larry R. Lyons, and David J. Gorney 

Space and Environment Technology Center, The Aerospace Corporation, El Segundo, California 


This is an investigation of stormtime particle transport that leads to formation of the ring 
current. Our method is to trace the guiding-center motion of representative ions (having selected 
first adiabatic invariants p) in response to model substorm-associated impulses in the convection 
electric field. We compare our simulation results qualitatively with existing analytically tractable 
idealizations of particle transport (direct convective access and radial diffusion) in order to assess 
the limits of validity of these approximations. For p i, 10 MeV/G (E H 110 keV at L « 3) the 
ion drift period on the final (ring-current) drift shell of interest (L ~ 3) exceeds the duration of 
the main phase of our model storm, and we find that the transport of ions to this drift shell is 
appropriately idealized as direct convective access, typically from open drift paths. Ion transport 
to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those 
portions of that drift path that lie outside the mean stormtime separatrix between closed and 
open drift trajectories. For p ~ 10-25 MeV/G (110 keV £ E < 280 keV at L « 3) the drift period 
at L « 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is 
transitional between direct convective access and transport that resembles radial diffusion. (This 
particle population is transitional between the ring current and radiation belt). For p £ 25 MeV/G 
(radiation-belt ions having E > 280 keV at L « 3) the ion drift period is considerably shorter than 
the main phase of a typical storm, and ions gain access to the ring-current region essentially via 
radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in 
this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of 
interest (L « 3), we have estimated (from both our forward and our time-reversed simulations) 
the time-integrated radial-diffusion coefficients D*jty for particles having selected values of p Z 15 
MeV/G. The results agree surprisingly well with the predictions (DlV) of quasilineax radial- 
diffusion theory, despite the rather brief duration 3 hr) of our model storm and despite the 
extreme variability (with frequency) of the spectral-density function that characterizes the applied 
electric field during our model storm. As expected, the values of deduced (respectively) 
from our forward and time-reversed simulations agree even better with each other and with D|JjP 
when the impulse amplitudes which characterize the individual substorms of our model storm are 
systematically reduced. 


1. Introduction 

This work is an outgrowth of our efforts to under- 
stand in detail the magnetospheric charged-particle 
transport that leads to stormtime ring-current formar 
tion. The ring current consists of geomagnetically 
trapped ions and electrons in the 10-200 keV energy 
range [e.g., Frank, 1967; Williams, 1981], and its global 
intensity is commonly measured by the geomagnetic 
index D 3t [e.g., Mayaud, 1980, pp. 115-129]. The 
quiet-time ring current is believed [e.g., Hamilton et al., 
1988] to contribute about 10-20 nT to -D s t, but this is 
largely offset by magnetopause currents during the ge- 
omagnetically quiet intervals which define the baseline 
(D s t = 0) for the index. Thus, the increase in \D s t\ to 
200 nT typically observed during the main phase of a 
major geomagnetic storm must correspond [cf. Dessler 
and Parker, 1959; Sckopke, 1966] to perhaps a ten- 
fold or twentyfold increase in the energy content of the 
trapped-particle population. 

Indeed, major increases in trapped-electron 
[ Williams and Smith, 1965; Pfitzer et al., 1966; Craven, 
1966; Soraas and Davis, 1968; Bostrom et al., 1970; 
Lyons and Williams, 1975] and in trapped-ion [Frank, 
1967; Frank and Owens, 1970; Smith and Hoffman, 
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1973; Williams and Lyons , 1974; Lyons and Williams, 
1976; Williams, 1981; Lui et al, 1987] fluxes are known 
to occur in connection with major geomagnetic storms. 
Such increases actually extend from L ~ 7 to as low 
as L ~ 2 and span particle energies from S 1 keV to 
several hundred keV. However, the main contribution 
to |D 3 i| comes from particles in the 10-200 keV energy 
range, since these contribute the overwhelming major- 
ity of the energy content of geomagnetically trapped 
particles. Trapped-particle enhancements at L £ 4 are 
commonly associated with isolated substorms, which 
can occur during otherwise quiet intervals, but the “in- 
jection” of ions and electrons to L S 4 corresponds to 
a geomagnetic storm, and particle injection to L 5 3 
corresponds to a major storm. Since particles injected 
to L < 4 are presumed to be the ones mainly respon- 
sible [Akasofu, 1963; Frank, 1967] for storm-associated 
increases in \D s t\, we must account for these flux in- 
creases in order to understand geomagnetic storms. 

Formation of the stormtime ring current presumably 
can be understood in terms of the transport of charged 
particles within the magnetosphere by storm-associated 
variations in the convection electric field. Early studies 
of charged-particle motion in the magnetosphere [e.g., 
Kavanagh et al., 1968; Chen, 1970] took into account 
the effects of a steady-state convection electric field. 
Chen [1970] described the adiabatic drift paths of equa- 
torially mirroring low-energy protons in a dipolar mag- 
netic field, with a radial corotation electric field and 
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a uniform steady-state convection electric field super- 
imposed. Later Ejiri [1978] traced particle trajectories 
under the influence of an increased steady-state con- 
vection electric field for arbitrary pitch angles. These 
models were able to account for substorm-associated 
enhancements of 1-100 keV particles observed near the 
plasmapause by Explorer 45. 

Although such steady-state descriptions are useful 
in explaining gross features of particle transport, tem- 
poral variations in the convection electric field must 
play an important role in both the injection and the 
trapping of particles during storms. Using a numerical 
model, Roederer and Hones [1974] discussed the char- 
acteristics of a time-varying electric field that could 
explain the frequent ATS observations [Deforest and 
Mcllwain, 1971] of plasma injection into the equato- 
rial magnetosphere. They found that low-energy par- 
ticles ( 5 10 keV) injected earthward from the tail by 
enhanced storm-associated electric fields can become 
trapped on closed drift paths when the enhanced elec- 
tric field decays away. Smith et al. [1979] illustrated 
this effect by tracing the motion of equatorial particles 
that had been injected from L = 10 on the night side 
in a dipolar field model. Later, Takahashi [1990], using 
the more realistic magnetic field model of Tsyganenko 
[1987], made a similar tracing of particle drifts. Injec- 
tion onto closed drift paths thus results directly from 
enhanced particle convection. 

The consequences of transporting low-energy ( 5 40 
keV) ring-current particles inward during a storm were 
discussed by Lyons and Williams [1980]. By mapping 
phase-space densities in accordance with Liouville’s the- 
orem, they showed that enhanced trapped-particle 
fluxes observed at L < 4 during two large storms could 
be accounted for by invoking a 1-3 Re inward radial 
displacement of the previously-trapped particle distri- 
bution under conservation of the first two adiabatic in- 
variants. The associated inward transport to a region 
of stronger magnetic field strength would result in the 
requisite particle energization. Assuming a duration 
~3 hr for the main phase of a storm, they found that 
this “displacement” would have required a mean en- 
hancement ~1 mV/m in the azimuthal component of 
the convection electric field across the nightside mag- 
netosphere. Particles having drift periods much longer 
than the main phase of the storm might then remain 
at a nearly constant local time (0) so that the en- 
hanced convection electric field could transport them 
uniformly inward. This model for particle transport 
into the ring-current region is called “direct convec- 
tive access” [Lyons and Williams, 1980]. It applies 
to low-energy particles ( $ 40 keV) and represents a 
limiting case for which the particle kinematics can be 
solved analytically. Particles with energies £ 40 keV 
would gradient-drift significantly in 0 during a 3-hr 
main phase and so would not satisfy the constant-0 
approximation. 

It is estimated, however, that 50-75% of the storm- 
time ring current energy is associated with particle en- 
ergies £40 keV [Williams, 1980]. Thus, an under- 
standing of the stormtime transport mechanism for 
higher-energy particles is at least equally important. 
Lyons and Schulz [1989] have recently treated the in- 
ward transport of £ 40 keV particles by using radial- 


diffusion theory as a first approximation. Using the 
random-impulse model of Cornwall [1968] to generate 
the fluctuation spectrum required in the resonant-par- 
ticle diffusion theory of Falthammar [1965], they found 
that storm-associated equatorial convection electric 
field fluctuations having root-mean-square magnitudes 
typical of those observed (~ 1.0-1. 5 mV/m) could ade- 
quately account for the increased trapped-particle 
fluxes typically observed at 2 ;$ L ;$ 4 during the main 
phase of a storm. Although radial-diffusion theory is an 
idealization of doubtful validity for ring-current particle 
energies S 100 keV, it represents a second limit which 
is analytically tractable. 

Direct convective access and radial diffusion are thus 
two limits that have been treated analytically and in- 
voked to account for the enhanced ring-current particle 
fluxes observed at 2 ;$ L S 4 during main phases of 
storms. However, much of the energy content of the 
ring current lies in an intermediate energy range (~40- 
100 keV) for which the assumptions invoked in order 
to arrive at these simplified analyses of the transport 
process are not valid. At such intermediate energies, 
the particle drift period can be comparable to the du- 
ration of the main phase of a storm, and so it is not 
clear how to describe the transport mechanism except 
via guiding-center simulations. 

We have undertaken such a guiding-center simulation 
to elucidate the process of stormtime ion transport and 
to assess the limits of validity of direct convective ac- 
cess and radial diffusion approaches. Our study entails 
certain innovations which facilitate the analysis of our 
results. 

1. We have synthesized an ensemble of model geo- 
magnetic storms, each consisting of almost randomly 
occurring exponentially decaying electrostatic impulses 
with quasi-random onset times. Since we express our 
electric field in analytical form, our approach enables 
us to make exhaustive diagnostic analyses on the con- 
sequences of each model storm in the ensemble. 

2. We have taken care to arrange our representa- 
tive test particles isochronally (i.e., at equal spacings 
in elapsed quiescent drift time rather than, for exam- 
ple, at equal spacings in magnetic longitude) around 
each drift shell of interest at the beginning of each sim- 
ulation. 

We have found that the stormtime access of ions hav- 
ing E % 110 keV to drift shells of interest in the con- 
text of ring-current formation is achieved typically by 
convection from the nightside neutral line along open 
drift trajectories. The minimum time required for this 
convective transport to occur is well approximated by 
generalizing the estimate of Lyons and Williams [1980] 
to the field model used in our study. The fraction of 
representative particles that gain access by direct con- 
vection to the drift shell of interest decreases with in- 
creasing ion energy, to the extent that particles having 
E £ 165 keV arrive almost exclusively from previously 
closed drift trajectories of either larger or smaller L 
value. For these particles we found that a transport 
description based on radial-diffusion theory is appro- 
priate. 

Indeed, for /x £ 15 MeV /G which corresponds to 
E £ 165 keV at L ~ 3 we have found interesting re- 
lationships upon comparing the radial-diffusion coeffi- 
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cients Dll deduced from our simulation results with 
the quasilinear predictions of Falthammar [1965] and 
Cornwall [1968]. In particular, we have found increas- 
ingly good agreement (even for an individual model 
storm having a 3-hr main-phase duration) between our 
simulation results and the quasilinear predictions of 
Falthammar [1965] upon decreasing the amplitudes of 
the electric field impulses in our ensemble of model 
storms. The quasilinear diffusion coefficient for an in- 
dividual storm is a remarkably variable function of the 
first adiabatic invariant /i, and we have found that 
Falthammar’s quasilinear theory accounts surprisingly 
well for the \i values at which our simulation yields espe- 
cially small values for Dll • Despite these encouraging 
numerical results, the use of quasilinear diffusion theory 
to describe radial transport is marginally justified from 
a fundamental standpoint at energies ( 1$ 165 keV) pri- 
marily responsible for the ring current. Thus, we have 
extended our simulations and our comparisons into the 
higher-energy radiation-belt range (ultimately to /L t = 
250 MeV/G, which corresponds to E = 2.75 MeV at 
L = 3). 

2. Field Model 

The magnetic field model used in this model is ob- 
tained [cf. Dungey, 1961; Hill and Rassbach, 1975] 
by adding a uniform southward field AB to the geo- 
magnetic dipole field. We invoke this simple field con- 
figuration because it enables us to make direct com- 
parisons between the simulated transport and previous 
analytical formulations. An advantage of our model 
over a purely dipolar field is the presence of a quasi- 
magnetopause at the boundary between closed and 
open field lines (see Figure 1). The equation of a field 
line in this model is 


L= L* 2L* 81* L-co 81* 2L* L= L* 



Fig, 1. An illustration of the magnetospheric magnetic field model 
used in this study. The model is symmetric about the sin 9 — 0 
axis and about the equatorial (z = 0) plane, which contains a cir- 
cular neutral line at r s» 6 on the magnetic shell L — L*. This de- 
fines a quasi-magnetopause which approaches an asymptotic dis- 
tance p* from the tail axis at large distances \z\ from the equatorial 
plane. 


1 + 0.5^0 (r/RE)csc 2 6 — const = L (1) 

where r is the geocentric distance, 6 is the magnetic 
colatitude, Re is the radius of the Earth, and b = 
1.5 L*Re = 12.82Re is the radius of the equatorial neu- 
tral line. This value of b, which is obtained by mapping 
the last closed field line (denoted L*) to a colatitude 
of 20° on the Earth, corresponds to |AB| = 14.474 nT 
and L* = 8.547. The limit b — » oo (L* — > oo) would 
correspond to a purely dipolar B field. In this study, 
we consider only equatorially mirroring particles. Our 
equatorial field intensity Bo is given by 

B 0 = - 14.474 nT (2) 

where y,E = 3.05 x 10 4 nT-J?|, is the geomagnetic dipole 
moment. Further details of this field model are given 
by Schulz [1991, pp. 98-110]. 

We assume that the total electric field E = — is 
derivable from the scalar potential 

in which the three separate terms correspond to coro- 
tation (Vjj = 90 kV), the Volland-Stern [Volland, 1973; 
Stem, 1973] model of quiescent convection (Vo = 50 
kV), and the time-dependent enhancement A V(t) as- 
sociated with the stormtime convection, respectively. 
The time-varying term in the potential is assumed to 
vary as L [cf. Nishida , 1966; Brice, 1967] rather than 
as L 2 because electric disturbances are expected to be 
less well shielded than steady-state convection by the 
inner magnetosphere. 

Whereas others have imposed time- varying cross- 
magnetospheric electric fields either estimated from Kp 
[e.g., Smith et al., 1979; Kistler et al., 1989] or deduced 
(insofar as possible) from direct observations [Riley and 
Wolf, 1992] during particular storms, we have instead 
synthesized an ensemble of model geomagnetic storms. 
In our model the storm-associated enhancement A V(t) 
in the cross-tail potential drop 

N 

&V(t) = X AVi ex P[(*i “ t )/ T } 0 ( t - ti) ( 4 ) 

i=l 

where B{t) is the unit step function (= 1 for t > 0; = 0 
for t < 0), as a superposition of almost randomly oc- 
curring impulses that rise sharply and decay exponen- 
tially with a “lifetime” r = 20 min [cf. Cornwall, 1968]. 
The impulses represent the constituent substorms of a 
storm. The potential drop A Vi associated with any 
impulse is chosen randomly from a Gaussian distribu- 
tion with a 200-kV mean and a 50-kV standard devia- 
tion. Negative impulses (such that AV) < 0) axe thus 
theoretically possible but extremely improbable, corre- 
sponding to only 0.003% of all impulses. We have cho- 
sen such a large mean for the AF; since our intention 
is to model a major (|D s t| ~ 200 nT) storm, such as 
those which Lyons and Williams [1980] analyzed. Since 
those storms had a main phase lasting ~3 hr, we assume 
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here that the N start times in (4) are randomly dis- 
tributed within a 3-hr time interval which corresponds 
to the main phase of our model storm. However, we im- 
pose a 10-min “dead time” (after each impulse onset) 
during which no subsequent impulse can start. This 
constraint imposes a realistic delay between the occur- 
rences of consecutive impulses and thus leads to a rea- 
sonable distribution of impulses. Without such a dead 
time it would be possible for the next impulse to start 
immediately after the previous one. We found (see Fig- 
ure 2a) that this could result in the buildup of unrealis- 
tically large cross-tail potentials, and so we rejected the 
model in which the f j in (4) were randomly distributed 
without dead time. The technical details of how we de- 
termine the start times of the impulses are given in the 
appendix. 

We have constructed 100 such random storms so that 
on average there are N = 9 impulses per storm or 3 sub- 
storms/hr. We have randomly chosen one model storm 
for a detailed case study. Figure 2b shows the varia- 
tion in cross-tail potential for this prototypical storm. 
The mean cross-tail potential drop over the time in- 
terval 1 1 < t < ti + 3 hr in this case is 230 kV, of 
which 50 kV corresponds to Vo and 180 kV to (A V(t)). 
This is a fairly typical value for (AF (f)) among the 100 
storms that we constructed since the ensemble-averaged 
(A V(t)} turned out to be 184 kV. 

3. Particle Dynamics 

Since we are simulating the guiding-center motion 
of nonrelativistic equatorially mirroring particles here, 
we treat the first two adiabatic invariants (fi ^ 0 and 
J = 0, respectively) as conserved quantities. In this 
case, the guiding-center Hamiltonian H is equal to the 
total energy (kinetic plus potential) of a particle. This 
is given by 

H = fiBo + ( 5 ) 

where q is the charge of the particle. It follows from 
(2)-(4) and from Hamilton’s equations of motion [e.g., 


Goldstein, 1980, p. 339] that the guiding-center motion 
of an equatorially mirroring particle subject to E x B 
and gradient- B drifts is described by the equations 


dL 

dt 


L 2 Re cos <(> 
2fi E 


v„(£)w<«)(£) 


( 6 ) 


and 


d<t> _ 0 

dt 6 qB 0 rS 






(7) 


where Ct is the angular velocity of the Earth and <p is the 
azimuthal coordinate (local time). We neglect particle 
loss processes such as charge exchange and Coulomb 
drag in this study. 

We solve the ordinary differential equations (6) and 
(7) simultaneously by using the Bulirsh-Stoer extrap- 
olation method with variable time step [e.g., Press et 
al., 1986, pp. 563-568] for specified initial conditions. 
First, by setting A F(i) = 0 in (6) and (7), we obtain 
steady-state adiabatic drift paths associated with a par- 
ticular value of the first adiabatic invariant p. We start 
12 representative particles at points equally spaced in 
time on a steady-state drift path. We then apply a 
storm-associated temporal variation AF(f) in the cross- 
tail potential drop and ! fun the simulation to determine 
the consequent stormtime particle transport. We can 
run the simulation either backward in time (to deter- 
mine where any representative particle must have been 
prior to the storm in order to reach the desired phase 
on its “final” drift shell) or forward in time (to follow 
the dispersal of initially codrifting particles among drift 
shells during and after the storm.) 



Fig. 2. (a) We have rejected the storm model in which the onsets of superimposed impulses in the cross- tail 
potential occur purely at random. Without any realistic delay (i.e., dead time) between $he onsets of consecutive 
impulses, it is possible to attain an unrealistically large cross- tail potential (e.g., ~ 650 kV at t « 1.6 hr in this 
instance). (6) We have adopted and applied a model storm in which the cross-tail potential is enhanced by 
a superposition of exponentially decaying impulses (decay time r = 20 min). These represent the constituent 
substorms of a storm and start at times which are distributed randomly over a 3-hr time interval (main phase) 
except that we impose, after the start of each impulse, a 10-min “dead time” during which no subsequent impulse 
can start. 
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The maximum time step in our simulations is re- 
strained so as to ensure good accuracy. We have tested 
our results against the conservation law [cf. Goldstein, 
1980, p. 348] 

dt dt ~ dR dt + d<j> dt ~ ° ^ ' 


on selected representative particles for time-dependent 
A V (t) as well as for Ah' (t) & 0 and have found (8) to 
be well satisfied throughout our simulations. 


4. Simulated Trajectories 

In this section we present results of the simulated 
transport of singly charged ions having various p values 
and vanishing second invariant (J = 0). The dashed 
outer circle (shown in whole or in part on each particle- 
trajectory plot) indicates the location of the neutral 
line, a circle of radius R = 12.82, which marks the 
boundary between open and closed magnetic field lines 
in our magnetic-field model (cf. Figure 1). The left- 
hand panels in Figures 3-4 show steady-state ion drift 
paths in the equatorial plane for selected n values. For 
fj, < 26.5 MeV/G there is an x-type stagnation point 
associated with the separatrix between open and closed 
drift trajectories. 

Figure 3a is for p, = 3 MeV/G, which corresponds to 
a kinetic energy of 33 keV at R = t/Re = 3. The ionic 
drift for this ji value has its stagnation point at R = 
4.4 on the dawn meridian. We label closed drift shells 
in terms of the dimensionless third adiabatic invariant 
defined by Roederer [1970, p. 107] as 


L = 


2 nfj,E 




( 9 ) 


where <5>s is the magnetic flux enclosed by that drift 
shell and L(</>) denotes the field-line label defined by 
(1) at longitude <f> on the drift shell. For example, the 
drift shell highlighted by the small solid circles in Figure 
3a bears the label L = 3.14. 

The right-hand panels in Figure 3-4 illustrate the 
stormtime access of particles to the drift path which 
crosses the dusk meridian at R — 3. This drift path has 
been chosen as representative of the particles that are 
primarily responsible [cf. Akasofu, 1963; Frank, 1967] 
for storm-associated increases in the ring-current index 
\D 3 t[. The steady-state drift period (T3 = 12.0 hr) for 
ions with n — 3 MeV/G on such a path is 4 times 
the duration of our model storm. “Starting” with 12 
representative particles (as indicated by the solid cir- 
cles) equally spaced in time on such a drift path, we 
have applied our model storm (see Figure 26) and have 
run the simulation backward in time so as to determine 
where the corresponding ions must have been located 
prior to the storm in order to have been transported 
to their final drift phases on the steady-state drift path 
in question. We followed open trajectories backward 
in time to the neutral line and closed trajectories back- 
ward in time for at least one quiescent drift period prior 
to the onset of the first impulse. To ensure that the 


last impulse has decayed by at least nine e-foldings, we 
initiate the time-reversed simulation at t = t\+ 3 hr 
+(m/12)r3 where T3 > 3 hr and m is an integer. We 
impose this last condition in a T3-dependent way sim- 
ply to make the representative ions’ phases correspond 
approximately to those shown in Figures 3-4. Since 
ions which differ in drift phase respond differently to 
the model storm, the right-hand panels of Figures 3-4 
illustrate (on an expanded scale) the range of possible 
responses. 

The small solid circles, spaced 30° apart in drift 
phase, in Figure 36 correspond to those in Figure 3a 
and thus mark the final steady-state drift path. Eight 
of the 12 representative particles could have been trans- 
ported inward from the night side to populate the final 
closed drift path. Only one of these 8 representative 
ions would have been located at L < L* = 8.547 at 
storm onset ( t = 0); its position at that time is indi- 
cated by the open circle. The other 7 of these 8 would 
have crossed the neutral line at times t > 0. Their ac- 
cess to the ring-current region is essentially convective 
[cf. Lyons and Williams, 1980]. Four of the 12 repre- 
sentative ions would have been transported outward to 
the drift shell of interest (L = 3.14) from closed drift 
shells of smaller L. 

Superposition of the mean stormtime separatrix 
(dashed-dotted curve) on the stormtime drift paths elu- 
cidates the condition for direct convective access to 
the drift shell of interest and further justifies the in- 
ference that this is the main mode of transport to the 
L = 3.14 drift shell for /a = 3 MeV/G. The mean storm- 
time separatrix is obtained by setting Vq = 50 kV and 
A V(t) = (A V(t)} = 180 kV in (3). It represents the 
boundary between open and closed adiabatic drift paths 
during the storm. Eleven of the 12 representative ions 
on the final steady-state drift path of interest lie out- 
side the region enclosed by the stormtime separatrix, 
and so we might expect approximately this fraction to 
have direct access to it during the storm. This is only 
a rough statistical prediction, however, and in fact our 
simulation showed that only 8 of the 12 representative 
particles gain access in this way. 

The case fj, = 10 MeV/G corresponds to an ion en- 
ergy of 110 keV at R — 3. Steady-state ion drift paths 
in the equatorial plane for this case are shown in Fig- 
ure 3c. The x-type stagnation point for /i = 10 MeV /G 
is located on the dawn meridian at R = 6.5, which is 
farther from the Earth than the stagnation point for 
H = 3 MeV/G. The drift shell that intersects the dusk 
meridian at R = 3 corresponds to L — 3.00 for n — 10 
MeV/G, and the quiescent drift period (T3 = 2.4 hr) 
for these ions has the same order of magnitude as the 
duration of the storm. 

Results of the time-reversed simulation for fi = 10 
MeV/G are shown on an expanded scale in Figure 3d 
For this higher value of n, the model storm has trans- 
ported fewer particles (only 6 of the 12) inward along 
open drift paths to the closed drift shell of interest in a 
manner which resembles direct convective access. The 
other half had already been on closed drift shells when 
the storm began, and 5 of these had actually started 
on drift shells interior to the final drift shell of interest. 
The final drift period is comparable to the duration of 
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the storm in this case, and the transport seems transi- 
tional between direct convective access (for half of the 
representative particles) and something that resembles 
radial diffusion for the other half. 

This point is further illustrated by superimposing the 
mean stormtime separatrix (dashed-dotted curve, Fig- 
ure 3d) on the final steady-state drift path. Because the 
stormtime separatrix lies slightly earthward of the final 
steady-state drift path of interest on the dusk side, par- 


ticles initially on open drift paths that pass near there 
can gain access during the storm to the region of closed 
steady-state drift paths. This argument should apply 
to only about half of the drift shell of interest and so 
the other half presumably could not have been popu- 
lated from open drift paths during this model storm. 
We cannot predict just which representative particles 
will gain access along open trajectories and which will 
not. Instead, about half of the L = 3.00 drift shell 


STEADY STATE 





Fig. 3. Quiet-time trajectories of singly charged ions having (o) /t = 3 and (c) 10 MeV/G are plotted in the left- 
hand panels. The outer dashed circle represents the neutral line at r = 6. Ions whose drift paths cross the dusk 
meridian at R = 3 have drift periods T3 = 12 hr and 2.4 hr for fi = 3 and 10 MeV/G, respectively, as noted. The 12 
representative ions’ “final” positions (corresponding to the beginning of our time-reversed simulation) on the final 
steady-state drift path of interest are marked by small solid circles (partially obscured in right-hand panels). The 
corresponding stormtime trajectories, as computed in our time-reversed simulation, are shown on an expanded scale 
in the right-hand panels (Figures 36 and 3d). The dashed-dotted curves represent mean “stormtime” separatrices. 
The small open circle indicates the “initial” position of the one representative ion that was on an open quiet-time 
drift trajectory but would have been inside the magnetosphere (i.e., at L < L*) at storm onset. The dashed portion 
of the corresponding trajectory applies to times before storm onset (t < ti). 
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Fig. 4. Quiet-time trajectories of singly charged ions having (a) p = 15 and( 6 ) 30 MeV/G are plotted in the left- 
hand panels. Ions which cross the dusk meridian at R = 3 have drift periods T 3 = 1.5 hr and 0.73 hour for p = 15 
and 30 MeV/G, respectively. The 12 representative ions’ “final” positions (corresponding to the beginning of our 
time-reversed simulations) on the final steady-state drift path of interest are marked by small solid circles (totally 
obscured in right-hand panels). Corresponding stormtime trajectories computed in our time-reversed simulation 
are shown on an expanded scale in the right-hand panels (Figures 45 and Ad). The dashed-dotted curves represent 
the mean stormtime separatrices for the respective p values. 


for p — 10 MeV/G should be populated by stormtime 
transport from open drift paths and about half from 
closed. Our simulation bears out this expectation. 

Ions having p = 15 MeV/G would have energies of 
167 keV at R = 3. Steady-state adiabatic drift paths 
for these (see upper left-hand panel of Figure 4) show 
a stagnation point at R = 7.4 on the dawn meridian. 
This is still farther from the Earth than for p = 10 
MeV/G. Again, the final drift shell of interest is that 
which intersects the dawn meridian at R — 3. The 


quiescent drift period T 3 for these ions is 90 min, which 
is only half the duration of the storm. 

The stormtime trajectories are plotted on an expand- 
ed scale in Figure 46. We find all the representative ions 
to have been transported from closed drift paths, some 
from smaller and some from larger L values. This leads 
to a spread among the initial L values. The stormtime 
transport for p = 15 MeV/G at L ~ 3 thus resembles 
radial diffusion although this was a time-reversed cal- 
culation. While Figure 46 clearly conveys the visual 
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impression of diffusive transport, the individual trajec- 
tories are not resolvable. 

Finally, we consider radiation belt-ions having p = 
30 MeV/G. These would have energies of 335 keV at 
R = 3. The quiescent drift period (73 = 44 min) of 
such ions on the shell that intersects the dusk meridian 
at R = 3 (and found to have L = 3.00) is less than 
1/4 of our model storm’s duration. A plot of the qui- 
escent adiabatic drift paths for p = 30 MeV/G (Figure 
4 c) shows no stagnation point earthward of the neu- 
tral line on the dawn meridian. When the storm is 
simulated, the time-reversed ion transport for p = 30 
MeV/G again resembles radial diffusion (see Figure Ad). 
About half of the representative particles would have 
reached the final drift path from higher L values and 
half from lower. More precisely, those particles that 
would have been on the dawn side at storm onset are 
the ones that have been transported inward since their 
initial gradient drift is toward the night side. Since 
the final steady-state drift paths lie well inside the re- 
gion bounded by the stormtime separatrix, there is little 
chance of stormtime access from an open drift path. 

We also ran simulations forward in time from t = 0 
for p Z 15 MeV /G so as to follow the dispersal of ions 
from a common initial drift path and in order to es- 
timate (see below) the corresponding radial-diffusion 
coefficient. Not surprisingly, the resulting stormtime 
transport was qualitatively similar to that shown in Fig- 
ure 4, although the particle populations followed in the 
time-forward simulations were implicitly different from 
those followed in the time-reversed simulations. This 
type of transport was typical for ions having p Z 13 
MeV/G, which corresponds to energies Z 145 keV at 
R = 3. 

Figures 3 and 4 thus illustrate a transition between 
modes of particle access to the stormtime ring cur- 
rent. Ring-current ions having p £ 3 MeV/G undergo 
direct convection from open (plasmasheet) drift paths 
to closed drift shells at R ~ 3, while radiation-belt ions 
having p Z 30 MeV /G respond to the same impulsive 
enhancement of the convection electric field in a man- 
ner that resembles radial diffusion among closed drift 
shells [cf. Lyons and Schulz, 1989]. The transition be- 
tween these two idealized modes of access occurs at p ~ 
5-15 MeV/G and thus (as might have been expected) 
for ring-current particles whose adiabatic drift periods 
are comparable to the duration of our model storm’s 
main phase. These transitional particles have energies 
~ 55-165 keV at R = 3 and are thus quite representa- 
tive [cf. Williams, 1981] of the ring current as a whole. 
Simplified analytical descriptions of the particle trans- 
port pertain most directly to these limiting cases (i.e., 
to p £ 3 MeV/G and to p £ 30 MeV/G, respectively), 
but application of the simplified analyses to interme- 
diate p values is similarly instructive for comparisons 
with our numerical results. Such comparisons are un- 
dertaken in the next two sections. 

Although our simulation results are presented in this 
paper as pertaining to singly charged ions, the results 
can be made applicable to multiply charged ions by 
replacing p in all our expressions with p/Z, where Z 
is the ionic charge state. For example, our results for 
singly charged ions having p = 3 MeV/G (shown in the 


upper panels of Figure 3) would apply equally well to 
doubly charged ions having p = 6 MeV/G. 

5. Direct Convective Access 
Lyons and Williams [1980] introduced a simple 
model for the direct convection of low-energy particles 
from one L value to another. In their model a particle 
is assumed to remain at fixed local time, while being 
transported radially by an uniformly enhanced azimu- 
thal electric field. Their model implicitly entailed a 
uniform convection electric field but can be generalized 
for application to our model. By treating cos <j> and AV 
as constant in time, we could integrate (6) to obtain the 
condition 


AVAt = 


heL* 
Re cos cj) 


_1 1 _ 2 ( 1 _ 

L f 2 Li 2 + k Ui 



■fc2 


In 


n + k/L f 

^ 1 + k/Li 


) 


( 10 ) 


for transport from Li to Lf (interpreted here as the 
initial and final field-line labels). The parameter k in 
(10) denotes (AV/Vo)L*. By taking the limit k — > 0, 
we recover from (10) the steady-state convection time 


2p g (L*) 2 
3VqRe cos (j> 




(11a) 


that would have been deduced directly from (6) by set- 
ting A V(t) = 0. By taking the limit k — ► 00 (Vo —* 0), 
we recover from (10) the storm-associated convection 
time 


A t = 


heL* 


Re cos <pAV \_L 2 f 




(116) 


obtained by Lyons and Williams [1980]. These last two 
expressions remind us that (10) should be applied only 
over time intervals during which cos cf> is reasonably con- 
stant. Application of (11) to longer time intervals would 
falsely predict eventual access to arbitrarily small Lf 
even for arbitrarily small Vo and/or AV, whereas coro- 
tation and/or gradient drift actually prevent this. 

Lyons and Williams [1980] had envisioned direct con- 
vective access from closed trajectories. However, our 
simulation results (see right-hand panels of Figure 3) 
indicate that many low-p particles reach R = 3 along 
open-drift trajectories and, in particular from L* during 
our model storm. We want to compare the transport 
times encountered in our stormtime simulations with 
predictions based on the simple model of direct convec- 
tion for these particles. Since AV varies with time in 
our model storm, we need to reinterpret AV(f)Af in 
(10) as the time integral of AV(<). We assume that 
<j> has been constant during most of the transport for 
low-p particles and thus replace <j> by <jf in (10), where 
4>* is the particle’s local time when it crossed the neu- 
tral line (L*) at time t — t*. We use the 180-kV value 
found from (AV(f)) during our model storm in order to 
treat the parameter k = (AV/Vo )L* as a constant on 
the right-hand side of (10). Thus, we obtain 



Chen et al.: Stormtime Transport of Ring Current and Radiation Belt Ions 


3843 


7(f) = AV{t)dt 
Jt* 

I1 E L* 1 1 _ 2(1 1 _\ 

Re cos <fr* Lj 2 7* 2 k\L* Lf) 




1 + k/ Lf 
1 + k/L* 


) 


( 12 ) 


as our generalization of (10). We plot the right-hand 
side of equation (12) against <f>* as the solid curve in Fig- 
ure 5. This curve represents our generalization of the 
convection-time estimate made by Lyons and Williams 
[1980], 

For comparison, we compute and plot 7(f) for each 
representative particle by integrating A V(t), as given 
by (4), over the time during which the particle was be- 
ing transported from the neutral line to the final drift 
path of interest. We found that 7 of the 12 represen- 
tative ions having /x = 3 MeV/G were transported to 
L = 3.14 in this manner. The corresponding values of 
7(f) are plotted as solid circles in Figure 5. Agreement 
with the direct-convective idealization is good for those 
ions which required the least 7(f), around 200-300 kV- 
hr, to get to Lf. Deviations from (12) are large for those 
ions which required larger values of 7(f), presumably 
because those ions have taken a more circuitous route 
to the final drift shell of interest. 



Fig. 5. The time-integral /(t*) = J 0 °° A V(t)dt required for con- 
vective transport of a particle from the neutral line (£*) to the 
final drift shell of interest, given by the right-hand side of (12), 
for fixed magnetic local time <fi = <j>*, is plotted as a solid curve. 
Solid and open circles represent values of I (t*) obtained from sim- 
ulated trajectories of representative ions having 3 MeV/G and 
10 MeV//G, respectively, under our model storm (cf. Figure 
2b). Solid and open squares represent values of /(t*) obtained 
from simulated trajectories of representative ions having 3 MeV/G 
and 10 MeV/G, respectively, under the idealized “boxcar” storm 
model specified by (21) with AV = 200 kV over a 3-hr main phase. 


As an additional test of (12) we make use of an earlier 
simulation we had run for ions having /x = 3 MeV/G in 
a storm model wherein AV(f) was represented as a box- 
car function of time (i.e., held constant at AV = 200 
kV over a 3-hr main phase). Values of the correspond- 
ing 7(f) = AVAt = (3 hr-f ) x 200 kV for those 7 of 
the 12 representative ions that were transported from 
the neutral line to the final drift path of interest during 
this idealized storm are plotted as solid boxes in Fig- 
ure 5. Agreement with (12) is best for those ions which 
had the shortest transport times (~1.2 hr) and thus the 
smallest values (~240 kV-hr) for 7(f). We have omit- 
ted plotting 7(f) in Figure 5 for the one representative 
ion having fi = 3 MeV/G that was on an open trajec- 
tory but located at L < L* at storm onset (see Figure 
36). This is because the right-hand side of (10) really 
applies with approximately constant k only to particles 
for which Li — L* , since k = 0 essentially by definition 
for t < 0. We expect that 7(f) for this particular repre- 
sentative ion would not agree well with (12) in any case, 
since its circuitous path leads to 7(f) ~ 600 kV-hr. 

We make similar comparisons in Figure 5 for ions 
having /x = 10 MeV/G. We plot values of 7(f) as open 
circles for those 6 of the 12 representative ions that were 
transported from the neutral line to the final drift path 
of interest during our applied model storm (see Figure 
26). Likewise, we plot values of 7(f) = (3 hr— f ) x 200 
kV as open boxes for those 5 of the 12 representative 
ions that were transported from the neutral line to the 
final drift path of interest during our boxcar idealization 
of a model storm. Again, agreement with (12) is good 
only for ions having short transport times ( t ~ 1.5 hr) 
in the boxcar idealization and relatively small values 
( ;$ 300 kV-hr) for 7(f) in general. However, a smaller 
fraction of the /x = 10 MeV /G ions than of the /x = 3 
MeV /G ions agree well with (12), probably because the 
effect of gradient-!? drift on the transport is greater for 
fj. = 10 MeV/G than for /x = 3 MeV/G. 

In view of the good agreement between values of 7(f) 
shown by small circles and by corresponding boxes in 
Figure 5, we infer that (12) will provide a good esti- 
mate for the minimum access time (for particles trans- 
ported most directly inward by mean stormtime AV) 
also in this situation. As in the case in Figure 5, how- 
ever, we expect to find that (12) will not provide a very 
good access-time estimate for particles that have taken 
a more circuitous route to the drift path of interest. 

6. Diffusion and Quasi-Diffusion Coefficients 

In contrast to the lower-/! ions, the transport of high- 
energy ions with fi £ 15 MeV /G (cf. Figure 4) resem- 
bles radial diffusion. In this section, we make detailed 
quantitative comparisons of radial diffusion coefficients 
obtained from our simulations with quasilinear theory 
[Falthammar, 1965] to show that the transport is in- 
deed well described as radial diffusion. For cases in 
which all 12 representative ions are found to have been 
transported from initially closed drift shells to the final 
closed drift shell which intersects the dusk meridian at 
R = 3, we have computed the L values Li (and Lf) 
of the respective initial (and final) drift shells in accor- 
dance with (9) and have constructed the quantity 

9 
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4=1 


(13) 

where T (= 3 hr) denotes the “duration” of the main 
phase of the model storm. The transport coefficient 
Dll defined by (13) corresponds to the familiar diffu- 
sion coefficient obtained in quasilinear theory 
[Falthammar, 1965], which applies in the limit 
maxjl-L^ 1 — LJ 1 \ — ► 0. The quantity DW“ is thus a 
measure of the variance among the initial third adia- 
batic invariants of particles situated on the final drift 
shell of interest. We also computed diffusion coeffi- 
cients D S £E from simulations run forward in time by 
interchanging the indices i and / in (13). The diffu- 
sion coefficients DffE obtained from simulations run 
forward and backward in time are not very different, 
80% of them being within a factor of 1.5 of each other, 
although they pertain implicitly to different sets of par- 
ticles. 

The resonant-particle formulation [Falthammar, 
1965; Schulz and Lanzerotti, 1974, pp. 81-95] of radial- 
diffusion theory leads to a diffusion coefficient of the 
form „ , 

L & R% 


Dll = 


4*4 


'-E 


(S> 


(14) 


where E(w/2w) is the spectral-density of the (quasi- 
uniform) equatorial electric field in the inner magneto- 
sphere and 1 * 3 /2tt is the particles’ drift frequency. Since 
the strength AE(t) of the fluctuating part of — is 
given, following (3), by 

Am --{- nk ^) Av ^ < 15 > 

in the equatorial plane at R -C L*, the spectral den- 
sity E(p 3 / 2 n) could alternatively be written as (1/2 L* 
Re) 2 V” (£ 13 / 2 - 11 ), where V{u/2-k) denotes the spectral- 
density function of A V(t), which is the storm-associat- 
ed enhancement in the cross-tail potential drop. The 
model for AV(t) specified by (4) leads to a spectral- 
density function of the form 


N N 


E(u/2n) — 


Tl + w 2 r 2 * 


Y Y AE i& E j cos[w(tj -ti)}, 


t=i j=i 


dead time constraint. A radial-diffusion model in which 
an infinite number of impulses (with an amplitude dis- 
tribution symmetric about zero) occur purely at ran- 
dom would lead to a statistical cancellation of the cross 
terms (j ^ i ) in (17). In this case, one would recover a 
radial-diffusion coefficient of the form found by Corn- 
wall [1968] and quoted by Lyons and Schulz [1989]: 

,_2 r 6 d2 N 

^ = WTixl(L*ni + n 2 3 T 2 ) fj (AKi)2 ‘ (18) 

The spectral density function E(w/2n) correspond- 
ing to our model storm (see Figure 26) is plotted as the 
solid curve in Figure 6. Its most notable feature is an 
extreme lack of smoothness [cf. Lanzerotti and Wolfe, 
1980]. The spectrum is so variable because the impulse 
onsets (i) in any individual storm modeled by (4) occur 
at specific times. The dashed curve in Figure 6 repre- 
sents the “smoothed” spectrum obtained from (16) by 
neglecting the cross terms (j ^ i) for the same model 
storm. We have introduced the smoothed spectrum in 
order to provide a standard of reference, equivalent to 
that introduced by Cornwall [1968], against which to 
compare our results for DJJE and Dj", (especially when 
these differ significantly from each otter). 

Since E(ui/2w) does not vary smoothly with lo, the 
quasilinear diffusion coefficient D^ L is not a smooth 
function of y. The solid curve in Figure 7 represents the 
variation of Df L with y for ions on drift shells which in- 
tersect the dusk meridian at R = 3. In order to plot the 
quasilinear diffusion coefficient as a continuous function 
of y, we used the empirical formulas 


73 : 


[1 + (0.17 //Lt) + (0.61 /m 2 )] x 10 4 
7.9434/u - 6.9418 * 


(19) 


and 


(16) 

where the summations extend over all N impulses (see 
the Appendix for details of the derivation). The corre- 
sponding diffusion coefficient given by quasilinear the- 
ory is 


, t 2 L 6 R% " A AVjAVj cos[n 3 (tj - tj)} 

LL 16Tm|(A*) 2 fa fa l + nb 2 

(17) 

The impulse-onset times ti and tj in our model storm 
are not quite randomly spaced in time because of our 


L = [0.335469 - (0.03/m) - 0.11/r 5/2 J - \ (20) 

which we inferred from our steady-state numerical re- 
sults by setting AV(t) = 0 in (6) and (7). Diffu- 
sion coefficients obtained from simulations run forward 
(solid circles) and backward (open circles) in time are 
also plotted in Figure 7 for comparison. Agreement of 
the diffusion coefficients obtained from the simulations 
with quasilinear theory is surprisingly good despite the 
strong variability of D q £ L with y. We find typically 
(i.e., in >80% of the available cases) that DTE and 
bf L agree within a factor of 4. Quasilinear theory even 
accounts for the y values (e.g., /x = 75 MeV /G and y 
— 80 MeV/G) at which the diffusion coefficients com- 
puted from the simulation are especially small. How- 
ever, there are still some discrepancies between 
and D^ l , typically at y values corresponding to sharp 
“peaks” or “valleys” in the electric spectral density. 
Cases in which the quasilinear diffusion coefficient un- 
derestimates D S EE (e.g., for y = 25 and 80 MeV/G) 
typically correspond to “valleys” in the spectral den- 
sity. Conversely, cases in which the quasilinear diffu- 
sion coefficient L>Y l overestimates DfE (e.g., for y = 
30, 60, 90, 100, 120, and 150 MeV /G) typically corre- 
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Fig. 6. The spectral density of the electric-field fluctuations in our 
model storm is plotted (solid curve) against normalized frequency 
(T = 3 hr). The dashed curve represents the “smoothed” spectral 
density computed as if the same impulses had been uncorrelated 
with respect to onset time (£*) and sign of amplitude (AV*). The 
dashed curve corresponds implicitly to the electric-impulse model 
considered by Cornwall [1968]. 

spond to “peaks” in the electric spectral density. For 
these cases, we find that the smoothed quasilinear dif- 
fusion coefficient (dashed curve in Figure 7) given 
by (18) generally agrees better with than does 
the unsmoothed quasilinear diffusion coefficient D g ^ L . 
In view of the extreme lack of smoothness in the spec- 
tral density this tendency of Djr to reduce the dis- 
crepancies between and D^ L led us to suspect 

that an appeal to resonance-broadening might improve 
upon our quasilinear results for Dj l [y The idea was that 


a spread in drift frequencies would accumulate among 
the 12 representative particles as a consequence of their 
spread in L. However, we found that such resonance- 
broadening only slightly smoothed the sharp minima 
and maxima in D q ^ L so that the resulting diffusion co- 
efficients D r £ L (n) still differed from Df£ by a factor 
£ 3 in about 20% of the cases [ Chen et al, 1992] in- 
stead of by a factor of £ 4 (as found here). 

Figures 8a and 86 show comparisons of the diffusion 
coefficients for model storms identical to that shown in 



Fig. 7. Solid and open circles represent values of radial-diffusion 
coefficient obtained via (13) from time-forward and time- 

reversed simulations, respectively. Quasilinear diffusion coefficient 
Df, (solid curve) is given by (17) and its “smoothed” counterpart 
D’l L (dashed curve) by (18) for model storm defined by Figure 
26. 




Fig. 8. Solid and open circles represent values of radial-diffusion coefficient obtained via (13) from time- 
forward and time-reversed simulations, respectively. Quasilinear diffusion coefficient Djf, (solid curve) is given by 
(17) as an implicit function of ft for model storms equivalent to that defined by Figure 2b except that the impulse 
amplitudes in A V(t) have been reduced by a factor of 2 (hence D^ L by a factor of 4) in Figure 8a and by a factor 
of 4 (hence D q £ L by a factor of 16) in Figure 86. 
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Figure 26, except that the impulse amplitudes A Vi were 
reduced by factors of 2 and 4, respectively. Agreement 
between Dy L and DYE improves considerably with the 
corresponding decreases in and (AV(t)), as we 
might have expected. The discrepancies between these 
estimates of the diffusion coefficients amount typically 
(i.e., in >80% of cases available for comparison) to a 
factor of 5 1.3 and S 1.2 when A V(t) is reduced by 
factors of 2 and 4, respectively. This agreement is re- 
markably good in view of the fact that the quasilinear 
diffusion coefficient D q £ L varies irregularly through 4 or- 
ders of magnitude for 15 MeV/G 5 ft ~ 250 MeV/G. 
Moreover, the agreement between diffusion coefficients 
obtained from simulations run forward (solid circles) 
and backward (open circles) in time improves notably 
as the impulse amplitudes in AV(f) are reduced. We 
find in >80% of the cases examined that the discrep- 
ancy between “time-forward” and “time-reversed” es- 
timates for DVE amount to factors ;$ 1.5, 5 1.3, and 
< 1.15 for (A v(t)) = 180, 90, and 45 kV, respectively: 
Improved agreement with decreasing (A V(t)) is not 
surprising since the quasilinear diffusion coefficient is 
formally invariant under time reversal. Quasilinear the- 
ory is supposed to approximate the simulated transport 
better and better as A V ( t ) — > 0, this being the limit in 
which the ion populations represented in simulations 
run forward and backward in time become less and less 
distinct. 

The surprisingly good agreement between D*£ L and 
DJ™ which we found for a storm modeled as a quasi- 
random sequence of impulsive enhancements in the con- 
vection electric field has led us to ask whether an even 
simpler storm model might lend itself to a successful 
quasi-linear description. As is noted above, we had suc- 
cessfully tested a boxcar representation, 

AV(i) = 0(t)9(T - t) AV, (21) 


with T = 3 hr, for the description of direct convective 
access for p ^ 13 MeV/G. We thought it might be in- 
structive also to compare the diffusion coefficients Df™ 
and D^r from the application of (21) with each other 
and with the diffusion coefficients plotted for p £ 15 
MeV /G in Figure 7. Thus, we have set A V = 200 kV 
in (21) for the purpose of making these comparisons. 
We might expect a much poorer agreement between 
the diffusion coefficients DYE and D q ^ L in this latter 
model for AF(t), but a formal substitution of (4) into 
(A2) via (16) nevertheless yields (see Appendix) 




( 22 ) 


Substitution of this spectral-density function into (14) 
yields a diffusion coefficient of the form 

* L e RU^V) 2 Tsm\nT/r 3 ) 

LL 16p|(L*)2 (xT/r 3 ) 2 ’ [ ) 

upon evaluation of the spectral density at the particle’s 
quiescent drift frequency fl^/2n s t^ - 1 . The quasi- 
linear diffusion coefficient is plotted as the solid curve 


in Figure 9. Because of the sinusoidal factor in the 
spectral-density function, the quasilinear diffusion coef- 
ficient vanishes for r 3 = T/(n + 1), with n = 0,1, 2, .... 
Thus, particles nominally executing an integral number 
of drift periods during the storm would experience no 
net diffusion, having been transported as much inward 
as outward during their n + 1 complete drift periods. 
For these cases the diffusion coefficient given by (23) 
will surely underestimate the actual transport coeffi- 
cient, since in reality there must be a spread in the 
stormtime drift frequencies of the representative parti- 
cles. The actual diffusion coefficient in such cases could 
possibly be estimated by averaging D q ^ L over an ap- 
propriate bandwidth of drift frequencies [cf. Dungey, 
1965; Schulz, 1975]. Here we smooth D q £* (perhaps ex- 
cessively) by replacing the factor sin 2 (7rT/r 3 ) with its 
mean value (1/2) so as to obtain 


.. L 6 R b (AV) 2 T3 
~ 32p|(L*) 2 7T 2 r 


for Ts <C \Ph\T « 800 min (this last inequality being 
well satisfied for the range of p values considered in 
this study). The dashed curve in Figure 9 represents 
D q £r- Since the boxcar function in A V(t) specified by 
(21) is symmetric in time about t = 1.5 hr, the diffusion 
coefficients DYE obtained from simulations rim forward 
or backward in time are equal. We compare diffusion 
coefficients obtained from the simulation (solid circles) 
with the quasilinear diffusion coefficient (solid curve) 
and with the smoothed quasilinear diffusion coefficient 
(dashed curve) in Figure 9. We find, quite remarkably, 
that D q ^ L actually provides a good estimate (within a 
factor of ~2) for DYE over the entire range of p values 
for which the uniform enhancement in the convection 



Fig. 9. Solid circles represent values of radial-diffusion coefficient 
D ] [JJ obtained via (13) from guiding-center simulations (either 
time-forward or time-reversed). Quasilinear diffusion coefficient 
dY (solid curve) is given by (23) and its “smoothed” counterpart 
DjEl (dashed curve) by (24) for idealized “boxcar” storm model 
denned by (21) with AV = 200 kV and T — 3 hr (main phase 
duration). 



Chen et al.: Stormtime Transport of Ring Current and Radiation Belt Ions 


3847 


electric field transports particles entirely among closed 
drift shells. ; .We find also (not surprisingly) that the 
values of D s £l . for ions at the larger p values (p > 90 
MeV /G) is generally greater during an impulsive model 
storm based on (4) than for a boxcar model based on 
(21). However, a comparison between Figures 7 and 9 
(note the one-decade scale offset) shows little system- 
atic variation between values of Df? deduced from (4) 
and from (21) for p £ 50 MeV /G. 

7. Summary and Conclusions 

We have used a dynamical guiding-center simulation 
model to study the stormtime ion transport which leads 
to formation of the ring current and diffusion in the 
radiation belts. We have traced representative ions’ 
guiding-center motion in response to model storm-asso- 
ciated impulses in the convection electric field for a 
range of ion p values. Although our simulation re- 
sults are presented in this paper as pertaining to singly 
charged ions, the results can be made applicable to mul- 
tiply charged ions by replacing p in all our expressions 
with p/Z , where Z is the ionic charge state. By com- 
paring our numerical results with previously formulated 
limiting idealizations of particle transport (convective 
access and radial diffusion), we have assessed the limits 
of validity of these approximations and have encoun- 
tered several interesting nuances. A summary of our 
results follows: 

, 1. For ions having drift periods that exceed the du- 
ration of the main phase of the storm (p < 10 MeV/G, 
hence E ;$ 110 keV at L ss 3), their inward transport 
to form the stormtime ring current is appropriately de- 
scribed as direct convective access [Lyons and Williams, 
1980]. 

2. For ions having drift periods comparable to the 
duration of the main phase of the storm (p ~ 10-25 
MeV/G, hence 110 keV <E< 280 keV at L « 3), 
there is a transition between direct convective access 
and transport that resembles radial diffusion. 

3. The mean “stormtime” separatrix constitutes an 
idealized boundary between closed and open particle 
trajectories for a given p. The location of the storm- 
time separatrix relative to the final (closed) adiabatic 
drift j>ath of interest has implications concerning the 
type of particle transport encountered. For example, 
any portion of the final drift shell of interest that lies 
outside the stormtime separatrix is accessible to par- 
ticles from previously open drift paths. Conversely, if 
the final drift shell of interest lies entirely earthward of 
the stormtime separatrix, then particles must be trans- 
ported from closed drift paths. 

4. For ions having drift periods considerably shorter 
than both T (the duration of the main phase of the 
storm), the transport is well described by radial diffu- 
sion. 

5. The electric spectral density of our model storm 
is extremely variable and leads to a highly irregular 
variation of radial-diffusion coefficients obtained 
from our simulations with p. The variation is similar to 
that found from observational data in spectral densities 
extracted by Lanzerotti and Wolfe [1980] and thus may 
be realistic. 


6. Despite their high variability with p, the diffu- 
sion coefficients Dfjf obtained from the simulations 
agree surprisingly well with the quasilinear estimates 
[cf. Falthammar, 1965] for Du at L « 3 for higher- 
p ions (p £ 15 MeV /G, hence E £ 165 kV at L « 3). 
Agreement between and the quasilinear diffusion 
coefficient ( D ¥,) improves when the impulse ampli- 
tudes in the convection electric field are decreased. 

An important implication of our study is that lower- 
energy (E ;$ 110 keV) ring-current ions at L « 3 are 
freshly injected there from open adiabatic trajectories, 
whereas the higher-energy (E S 150 keV) ring-current 
population consists of a mixture of freshly injected and 
previously trapped ions (both having been transported 
quasi-diffusively to the drift shell of interest). Composi- 
tion measurements on the ring-current particle popula- 
tion region before and after the main phase of a geomag- 
netic storm might provide a test of this interpretation. 

APPENDIX: 

Construction and Analysis of Model Storm 

In order to determine a reasonable set of f j for use 
in (4) to construct our model storm, we have gener- 
ated 1800 random real numbers between 0 and 100 and 
have sorted these into 100 bins by noting the integer 
part of each random number. Although our proce- 
dure strictly (rather than just probabilistically, as it 
should have) constrains the total to 1800 chosen ran- 
dom numbers per 100 bins and thus the mean to 18 
chosen random numbers per bin, we have found that 
the distribution of random numbers among bins is well 
approximated by a Poisson distribution with a mean 
of 18 per bin. In other words, the number of bins con- 
taining N chosen random numbers is well approximated 
by 100(N n /N\) exp(— N) with N =18. A chi-square 
test of our distribution of random numbers against a 
Poisson distribution with N = 18 yields x 2 — 9.50 for 
14 degrees of freedom. The probability P that any al- 
ternative distribution of random numbers would have 
yielded a smaller deviation from the expected bin val- 
ues for 14 degrees of freedom is P(x 2 < 9-50) = 0.20 
[see Hogg and Craig, 1959, p. 237]. Since this P lies 
well between 0.1 and 0.9, the x 2 test can be interpreted 
as having confirmed a Poisson distribution with N = 18 
[cf. Evans, 1955, p. 775]. 

Each bin represents a 3-hr interval from the ensem- 
ble of model storms. Our intent is that, on average, 9 
actual impulses (identified by the tf) should occur dur- 
ing the “typical” main phase, in which case half of the 
“typical” 3-hr interval would be disqualified by “dead 
time” from containing a tp We multiply the fractional 
part of each chosen random number by 180 min and 
identify the earliest time thus generated in each storm 
of the ensemble with i = 1 in (4). Since we impose a 
10-min “dead time” after each t, in (4), we reject sub- 
sequent random times that exceed any such <j by less 
than 10 min. In other words, i = 2 in (4) corresponds 
to the next random number that exceeds the smallest 
random number in the same bin by at least 1/18, and 
so on. Thus, the distribution of actual impulse start 
times ti is significantly distorted relative to a Poisson 
distribution with N = 9 [cf. Muller, 1973, 1974], even 
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though there are 9 impulses in the “typical” 3-hr main 
phase or 3 substorms/hr on average. For example, a 
Poisson distribution with N = 9 would have yielded 38 
of 100 storms with 9± 1 impulses, whereas we find that 
77 of our 100 randomly generated main phases contain 
9 ± 1 “actual” impulses. Indeed, our distribution cor- 
responds more closely to a Muller distribution, based 
on a Poisson distribution with N = 18 but distorted by 
a 10-min “dead time” [cf. Muller , 1973, 1974], accord- 
ing to which we should have expected 67 of 100 storms 
to contain 9 ± 1 “actual” impulses. A chi-square test 
of our distribution against a Muller distribution yields 
X 2 = 4.13 for 5 degrees of freedom. The probability 
P that any alternative distribution of random numbers 
would have yielded a smaller value of x 2 for 5 degrees 
of freedom is P(x 2 < 4.13) = 0.47, and so the x 2 test 
confirms [cf. Evans, 1955, p. 775] that our ensemble 
of model storms satisfies the appropriate Muller dis- 
tribution. Since our dead time constraint disqualifies 
randomly chosen candidate start times that follow too 
closely upon actual impulses, we find that the spacing 
between actual impulses tends to be more nearly even 
than random. 

We have constructed our ensemble of model storms to 
have (on average) 9 impulses with 20-min decay times 
during a 3-hr main phase, we should have expected the 
time integral of A V(t) from t\ to oo to be ~600 kV- 
hr. Indeed, the time integral of AV(t) from tj to oo 
for our prototypical storm is 594 kV-hr. We have noted 
above that (AV(f)) for our model storm is 180 kV. The 
ensemble-averaged (AV (i)) among the 100 storms that 
we constructed turned out to be 184 kV. This means 
that the time integral of AV(t) from t\ to ti + 3 hr 
would have been 540 kV-hr and 552 kV-hr for our model 
storm and the ensemble, respectively, but this 3-hr in- 
terval of integration excludes a significant portion of 
the final impulse. 

Quasilinear radial-diffusion theory [Falthammar, 
1965] requires us to obtain the spectral density function 
of the electric field associated with the realization of our 
model storm that is illustrated in Figure 26. Blackman 
and Tukey [1958, p. 7] define the spectral density func- 
tion E(ui/2tt) as follows: 


E(ul2ix) = lim 
T— > oo 


f\i: 


■T/2 1 2 

A E(t) exp (—iu>t)dt 

T/2 


(Al) 


for — oo < to < oo. We need to modify this definition 
for use in (4) because we regard the diffusion coefficient 
as a function of time and presume it to be negligible for 
t < 0 and for t > 3 hr. However, we would not wish to 
truncate the integral in (Al) at time t = 3 hr, since the 
contribution of the last few impulses to particle trans- 
port is not negligible until somewhat later. We truncate 
our simulations of particle transport during our model 
storm around t = 6 hr or later, but that choice is quite 
arbitrary. 

There would be no dilemma concerning the appro- 
priate definition of E(w/ 2n) if the dynamical process 
being modeled were statistically stationary [cf. Black- 
man and Tukey , 1958, pp. 4-7], but in that case we 
would not be modeling a storm. We adopt the defini- 
tion 


1 f°° 

E(fji/ 2 tt) = — / AE(t)exp(-itvt)dt 

* J — OO 


(A2) 


with T (= 3 hr) being the span of time over which im- 
pulse onsets can occur, as the compromise appropriate 
to the present scenario. 

The quantity actually simulated in the analysis pre- 
scribed by (13) is not the instantaneous Du, but rather 
its time integral, which corresponds to the ensemble 
of ion transitions between initial and final L values. 
We have identified this time integral with TD^™. Ac- 
cording to the classical analogue of Heisenberg’s uncer- 
tainty principle, it would not be meaningful to specify 
E(u/2n) as a continuous function of both frequency 
and time. In effect, however, a comparison of our re- 
sults for DYi with quasilinear theory would amount 
to a comparison of TDf£ with the time integral of 
(L 6 f|£/4^| ; )J&(fi3/27r), which we identify via (14) with 
TD\ h . In other words, our division by T in (14) and 
(A2) amounts to a convention that facilitates compar- 
isons with data which are best interpreted in terms of 
a time-dependent radial-diffusion coefficient [cf. Lanze- 
rotti et al., 1971; 1978; Lyons and Schulz , 1989], where- 
as (A2) formally enables Dfr to be calculated via (14) 
as a continuous function of drift frequency. For exam- 
ple, the electric-potential model A V(t) specified by (14) 
leads to a spectral-density function of the form 


hi yoo 

E(u/2n) = ~ X / A Ei exp[r _1 (tj - 1) + iojt\dt 
1 «= i J *i 


1 T 2 

Tl+w¥ 


N N 


X X AE i AE j COS [u{tj - ti)), 

i=l j = 1 


(A3) 


where the summations extend over all N impulses. The 
application of (A3) to the model storm illustrated in 
Figure 26 leads to the spectral-density function shown 
in Figure 6 and to the radial-diffusion coefficient plotted 
in Figure 7 as a continuous function of \i. 
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